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We report the discovery of a highly coherent oscillation in a type-I X-ray burst ob- 
■ served from 4U 1916—053 by the Rossi X-ray Timing Explorer (RXTE). The oscillation 

(N 
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it increased in frequency over the following 4 seconds of the burst decay to a maximum 
of ~272 Hz. The total measured drift of 3.58 ± 0.41 Hz (la) represents the largest 
fractional change in frequency (1.32 ± 0.15 %) yet observed in any burst oscillation. 



If the asymptotic frequency of the oscillation is interpreted in terms of a decoupled 



■ surface burning layer, the implied neutron star spin period is around 3.7 ms. However, 

the expansion of the burning layer required to explain the frequency drift during the 
burst is around 80 m, substantially larger than expected theoretically (assuming rigid 
rotation). The oscillation was not present in the persistent emission before the burst, 
nor in the initial rise. When detected its amplitude was 6-12% (RMS) with a roughly 
sinusoidal profile. The burst containing the oscillation showed no evidence for photo- 
spheric radius expansion, while at least 5 of the other 9 bursts observed from the source 
by RXTE during 1996 and 1998 did. No comparable oscillations were detected in the 
other bursts. A pair of kilohertz quasi-periodic oscillations (QPOs) has been previously 
reported from this source with a mean separation of 348 ± 12 Hz. 4U 1916—053 is the 
first example of a source where the burst oscillation frequency is significantly smaller 
than the frequency separation of the kHz QPOs. 

Subject headings: accretion, accretion disks — stars: individual (4U 1916—053) - 
stars: neutron — X-rays: bursts 



1. Introduction 



Type I X-ray bursts are signatures of unstable nuclear burning on the surface of accreting 
neutron stars (see Lewin, van Paradijs, & Taam 1995 and Bildsten 1998 for recent reviews). Strong 
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evidence for surface brightness anisotropics during the bursts comes from observations of highly- 
coherent 'burst oscillations' in 7 sources to date (Strohmayer et al. 1996; see also van der Klis 2000 
and references therein). These oscillations have been suggested to result from initially localized 
nuclear burning, which spreads over the surface of the neutron star during the early stages of 
the burst (Strohmayer et al. 1996). As the burst evolves, the oscillations typically increase in 
frequency by a few Hz; in X1658— 298 the increase is « 5 Hz, the greatest increase seen so far 
(Wijnands et al. 2000). This frequency drift may occur as a consequence of the burning layer 
becoming decoupled from the star (Strohmayer et al. 1997; Cumming & Bildsten 2000). In several 
sources the oscillation approaches an asymptotic frequency that may be the spin frequency f sp i n 
of the neutron star. Although many observational characteristics of burst oscillations seem to be 
consistent with this picture, a growing body of observations points to a substantially more complex 
underlying mechanism (e.g. Miller 1999; Strohmayer et al. 2000). In this Letter we report the 
discovery of a new burst oscillation in 4U 1916—053. 

The low-mass X-ray binary (LMXB) 4U 1916-053 (/ = 31.4°, b = -8.5°) was discovered by 
the Uhuru satellite in 1977 (Forman et al. 1978). EXOSAT observations revealed irregular X-ray 
dipping behavior with a period of ~ 50 min (Walter et al. 1982; White & Swank 1982), which 
optical observations of the my = 21 companion confirmed was approximately the orbital period 
(Grindlay et al. 1988). Analysis of RXTE observations during 1996 shows that the source is a 
member of the 'atoll' class (Bloser et al. 2000; Boirin et al. 2000). Both high (> 100 Hz) and 
low frequency QPOs were detected in the power spectra of the source as measured by RXTE. In 
particular, a pair of high-frequency QPOs ('kilohertz QPOs') with a mean separation of 348±12 Hz 
was simultaneously detected on five occasions (Boirin et al. 2000). On four of these occasions, the 
observed separation was consistent with the mean value, but on one occasion was only 290 ± 5 Hz. 
X-ray bursts exhibiting photospheric radius expansion resulting from super-Eddington luminosities 
indicate a source distance of 8.4-10.8 kpc (Smale et al. 1988). The X-ray bursts observed by RXTE 
during 1996 were searched for high-frequency oscillations, but none were found to a 3a limit of 3% 
RMS amplitude (Boirin et al. 2000). Here we present the results of a more comprehensive timing 
study of X-ray bursts in 4U 1916—053 observed by RXTE, including both the previously analyzed 
1996 data and more recent observations made during 1998. 



2. Observations and Analysis 

4U 1916—053 was observed using the RXTE Proportional Counter Array (PCA; Jahoda et al. 
1996) throughout 1996 and during 1998 June-August. These observations comprise 40 separate 
pointings with a total exposure time of 323 ks. The data were screened to exclude earth occulta- 
tions and intervals of unstable pointing, with the maximum allowed offset 0.01°. A search of the 
Standardl mode data (2-60 keV, 0.125 s time resolution, no energy resolution) revealed four type 
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I bursts during 1996 and six during 1998 1 (see Table 1). It is possible that one or more of the 
observed bursts originated from a source other than 4U 1916—053 in the PCA field (FWHM = 1°). 
The nearest (~ 1°) known source from the Simbad catalog (http://simbad.u-strasbg.fr/Simbad) 
is RXS 192242.1-051559, from which no bursts have been reported. We consider the likelihood of 
an unknown bursting source within the PCA field of view at this relatively high Galactic latitude 
remote. As already noted by Chou et al. (2000), the 1996 bursts are clustered in binary phase 
near the X-ray dip. Interestingly, the 1998 bursts do not show the same evidence for clustering, 
although the Chou et al. ephemeris still correctly predicts the X-ray dip times in the 1998 data. 

The burst profiles typically exhibit a broad and somewhat flat-topped peak, often followed 
by one or more narrow dips during the decay. Uninterrupted PCA data from high-time-resolution 
modes (typically the E_125us_64M_0_ls mode with 122 /is time resolution and 64 channel energy 
resolution) were available for 8 of the 10 bursts observed. The spectral characteristics of these bursts 
were measured by fitting absorbed blackbody spectra in the 2.5-40 keV range to 0.25 s segments 
of data covering the burst. Spectral fitting was undertaken with xspec version 11 (Arnaud 1996). 
A spectrum extracted from a short interval of data immediately before the burst was used as the 
background. Five of the bursts showed an increase in the fitted blackbody radius coupled with a 
decrease in the effective temperature during the first 10 s. This strongly suggests that photospheric 
radius expansion occurred during these bursts. No evidence for radius expansion was found for the 
bursts on 1998 July 24, August 1, and August 10. 

Sub-intervals of the high-time-resolution data covering each of the bursts were searched for 
oscillations in the 0-4000 Hz range. We computed 8x oversampled Fourier power spectra for 0.5 s 
lengths of data spaced at 0.25 s intervals and searched for statistically significant pulsed signals. A 
strong, highly significant (chance probability for 2000 independent frequencies 5 x 10~ 6 , equivalent 
to 4.6<t) oscillation at ^burst ~ 270 Hz was detected 0.75 s after the burst rise during the 1998 August 
1 burst (Figure 1). Weaker oscillations were measured at a gradually increasing frequency over the 
next ~ 5 s, except for the interval between 2 and 3 s following the burst rise when no significant 
oscillation was detected in the 220-320 Hz range (Figure 2a). When the oscillation reappeared 3 s 
after the burst rise, we initially detected two peaks of approximately equal strength separated by 
~ 3 Hz, although clearly only one was consistent with the continuing evolution of the initial signal. 
We note that the presence of two closely separated peaks has also been observed in 4U 1636—536 
(Miller 2000). The oscillation was not detected in the persistent emission prior to the burst, nor 
during the burst rise itself. The frequency and RMS amplitude for the oscillation corresponding to 
each peak was measured from the oversampled Fourier power spectra (e.g. Middleditch & Nelson 
1976). The difference between the initial and final frequencies at which the oscillation was detected 
(268.45 ±0.34 and 272.03 ±0.24 Hz respectively) was 3.58 ±0.41 Hz (la error; Figure 2b). The RMS 



1 An apparent weak burst on 1996 March 13 was on closer examination present in only one of the five proportional 
counter units (PCUs). We attribute this to a detector breakdown in PCU 3, the earliest example of which was 
previously thought to occur several days later. 
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amplitude for the full PC A energy range (2-60 keV) peaked at &12% around 1 s following the burst 
rise, and then decreased within 0.5 s to unmeasurable levels. When the oscillation returned ~ 3 s 
following the burst rise, the measured amplitude was in the 6-11% range. The RMS amplitudes 
were also calculated in the 2-8 and 8-20 keV energy ranges separately. The amplitude was generally 
significantly greater in the 8-20 keV energy band, at 6-12%. However, at times the amplitude in 
the 2-8 keV band reached that of the higher energy band, suggesting spectral variations in the 
pulsed emission. We note that when two significant peaks close to 270 Hz were measured in the 
Fourier spectrum (at around 2.5 s following the burst onset) only the higher frequency oscillation 
was present to a significant level in the 8-20 keV energy band. The pulse profile was generally 
consistent with a sinusoid. 

No significant persistent oscillations were observed in the other 9 bursts detected from the 
source in 1996 and 1998. A more sensitive search was undertaken on those bursts by averaging 
power spectra over 1 s segments of data within intervals of varying length covering the burst peak, 
but this too resulted in no detections. In several of the bursts (1996 August 16 in particular), 
peaks representing detections at > 90% confidence around 256 or 270 Hz were found in single 
power spectra as early as 2 s prior to the burst or coincident with the burst rise. These peaks did 
not persist for more than 0.5 s, nor were they as significant in power spectra of longer stretches 
of data. The most significant peak (3a equivalent for 4000 independent frequencies) was found in 
both the 0.5 and 1 s power spectra around 4 s before the 1996 August 16 burst, with frequency 
278.01±0.07 Hz (la) and estimated RMS amplitude 24%. We include these results for completeness; 
without repeated detections or greater significance their authenticity is questionable. 



3. Discussion 

We have conclusively detected a highly coherent ~270 Hz oscillation during a type I X-ray 
burst from 4U 1916—053, bringing the total number of known burst oscillation sources to eight. 
If we interpret the frequency evolution of the oscillation during the burst in terms of a decoupled 
burning layer, then the maximum observed frequency implies a neutron star spin period of ~ 3.7 ms. 
However, the large change in frequency during this burst (3.58 ± 0.41 Hz or 1.32 ± 0.15 %) would 
require an « 80 m expansion of the burning layer (assuming rigid rotation). The maximum peak 
flux of all the bursts observed is a lower limit to the Eddington flux -FecM for the source, in which 
case the peak flux for the 1998 August 1 burst (in which the burst oscillation was observed) is 
at most 0.5.FEdd- The maximum expansion predicted for a mixed H/He burst at this flux level 
is only 20 m (Cumming & Bildsten 2000), different from our measurement at the 7a level. That 
the binary period of the source is shorter than 80 min indicates that the mass donor must be very 
hydrogen-poor (Nelson et al. 1986), suggesting a predominantly He burning layer. Interestingly, 
Cumming & Bildsten (2000) predict an even smaller expansion in that case. The limited sample 
available suggests that burst oscillations in 4U 1916—053 are more likely to be present in atypically 
weak bursts which do not exhibit photospheric radius expansion (see Table 1). The other two 
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sources with ~ 300 Hz burst oscillations (4U 1728-34 and 4U 1702-429) also show this tendency, 
while the sources with higher frequency oscillations behave in a distinctly different manner (Muno 
et al. 2001, in preparation). 

In the four burst sources from which both a kHz QPO pair and a burst oscillation were 
previously known, the burst oscillation frequency ^ burst is comparable to (or else nearly twice) the 
peak separation of the kHz QPOs Av (see Table 2). Precise measurements of Av have been made 
for two of these sources, and these detailed measurements show clearly that 

^burst>"Az/ with n = 1 or 2. (1) 

(Mendez et al. 1998; Mendez & van der Klis 1999). For the other two sources, one can only conclude 
that Au and i^burst ar e equal within their uncertainties. The relationship between z^burst and Av 
given in (1) has been interpreted in terms of a beat frequency model for the kHz QPOs (Strohmayer 
et al. 1996; Miller, Lamb, & Psaltis 1998). In the original sonic-point beat frequency model of Miller 
et al. (1998), the upper kHz QPO appears at the Keplerian frequency of the accretion flow's sonic 
point, and the lower kHz QPO appears at the beat between this frequency and v s . In this simple 
picture, we expect Av = z^pm- More recent calculations by Lamb & Miller (2000) have shown 
that the accreting material's inward drift decreases the upper kHz QPO frequency below the sonic 
point's Keplerian frequency and increases the lower kHz QPO frequency above the beat frequency, 
resulting in Av < f sp i n - This is consistent with the observed relation (1) if we have i^burst ~ n ^spin 
with n = 1 or 2, as in the decoupled burning layer model for the burst oscillations. 

Our study of 4U 1916—053 has revealed the first example of a source where i^burst < Av 
to high significance (> 7a for the range of values of Av reported in Boirin et al. 2000). It is 
unclear how to understand this result in terms of the sonic point model, which would require the 
frequency of the upper kHz QPO to be higher than the Keplerian orbital frequency in order to 
match our observation. The numerical simulations presented by Lamb & Miller (2000) do not show 
this trend, but such a possibility cannot be excluded a priori. A different explanation for the kHz 
QPOs involving relativistic precession of the accretion disk has also been proposed. Here, the upper 
kHz QPO appears at the Keplerian frequency at some characteristic disk radius, and the lower kHz 
QPO occurs at the periastron precession frequency at the same radius (Stella & Vietri 1998; Psaltis 
& Norman 2000). In this picture, Av is approximately equal to the radial epicyclic frequency and 
is not inherently related to f S piru so that the striking observed similarity between z^burst and Av 
(or 2Az^) is not explicitly addressed. An additional model is therefore required to understand the 
burst oscillations in conjunction with the relativistic precession models for the kHz QPOs. 

The requirements for such a model have been discussed, although none has yet been developed 
quantitatively. If we assume that z^burst ~ ^spim then it must be that v sp i n is driven towards 
the maximum epicyclic frequency (Stella 1999) or else that the kHz QPOs only reach detectable 
amplitudes when Av ps v sp i n or Av v spm /2 due to a resonance effect (Psaltis 2000) Alternatively, 
z^burst might not be related to v sp i n at all, but might instead be related to an accretion disk mode 
that occurs near the maximum epicyclic frequency at the innermost stable orbit (Titarchuk et al. 
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1998; Psaltis & Norman 2000), such as the (/-modes predicted in the case of accreting black holes 
(see e.g. Kato 1990; Nowak & Wagoner 1991). Unfortunately, the absence of a detailed link between 
the burst oscillations and the relativistic precession model makes it impossible to constrain this 
model using our results. 

We thank Ed Morgan and Dimitrios Psaltis for useful discussions, and Josh Grindlay for 
making a source ephemeris available prior to publication. This work was supported in part by the 
NASA Long Term Space Astrophysics program under grant NAG 5-9184. 
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Fig. 1. — Leahy-normalized power spectrum from RXTE observations of 4U 1916—053 on 1998 
August 1. The 8-times oversampled power spectrum is calculated from 0.5 s segments of data, 
with frequency resolution 0.25 Hz. a) Power spectrum from 0.5 s segment beginning 0.75 s after 
the start of the burst. The largest peak is at a significance level of 99.99% and indicates a burst 
oscillation at 269.4 Hz. Expected positions of the n = 2,3,4 harmonics and the first subharmonic 
are indicated by the arrows, b) Power spectrum for the 0.5 s segment beginning 2.5 s after the start 
of the burst. The less significant marked peak is close to the (single) peaks observed in subsequent 
power spectra. 
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Fig. 2. — Time evolution of the oscillation in the 1998 August 1 burst from 4U 1916—053. (a) Dy- 
namic power spectrum of the burst oscillation (contours) overplotted on the burst intensity profile 
(solid line). The contour levels start at a Leahy- normalized power corresponding to 90%-confidence 
and increase in steps of 5 in power. The highest contour level corresponds to a confidence level of 
> 99.99%. (b) Time evolution of the oscillation frequency and RMS amplitudes. The top panel 
plots ^burst calculated from the oversampled Fourier transforms. The three panels below give the 
RMS amplitude in the full PCA range (2-60 keV) and the sub-ranges 2-8 and 8-20 keV respectively. 
Error bars represent the la uncertainties. 
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Table 1. X-Ray Bursts from 4U 1916-053 





Start time of burst 


Peak flux a 


Radius 


Oscillation 


Binary 


Observation ID 


(UTC) 


(10~ 8 erg cm -2 s _1 ) 


expansion? 


frequency (Hz) 


phase b 


10109-01-04-00 


1996 May 5 22:18:08 


2.17 ±0.26 


? 




0.02 


10109-01-05-00 


1996 Jun 1 18:28:44 


2.90 ±0.23 


? 




0.92 


10109-01-07-00 


1996 Aug 16 12:27:37 


3.06 ±0.17 


Y 




0.04 


10109-01-09-00 


1996 Oct 29 06:57:51 


2.72 ±0.14 


Y 




0.02 


30066-01-02-07 


1998 Jul 23 05:30:52 


3.26 ±0.18 


Y 




0.12 


30066-01-02-08 


1998 Jul 23 11:44:37 


3.34 ±0.19 


Y 




0.59 


30066-01-03-00 


1998 Jul 24 17:03:15 


2.23 ± 0.24 


N 




0.76 


30066-01-03-02 


1998 Jul 26 17:05:16 


3.16 ±0.17 


Y 




0.38 


30066-01-03-03 


1998 Aug 1 18:23:45 


1.62 ±0.13 


N 


269-272 


0.71 


30066-01-03-04 


1998 Aug 10 11:40:53 


3.03 ±0.15 


N 




0.79 



a Bolometric flux calculated according to the blackbody spectral fit parameters 

b Phase with respect to 50 min X-ray dip period, with dip phase at 0.0. From ephemeris of Chou et al. 
(2000). 
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Table 2. Burst Oscillation Frequencies and kHz QPO Peak Separations 



Source 


V burst (Hz) 


Ai/ (Hz) 


Refs. 


4U 1636-53 


291, 582 


251(4) 


1, 2,3 


4U 1702-43 


330 


315(ll)-344(7) 


4, 5 


4U 1728-34 


363 


279(12)-349(2) 


6, 7 


KS 1731-260 


524 


260(10) 


8, 9 


4U 1916-053 


270 


290(5)-350(25) 


10, 11 



References. - (1) Strohmayer et al. 1998; (2) Miller 
1999; (3) Mendez et al. 1998; (4) Strohmayer & Mark- 
wardt 1999; (5) Markwardt et al. 1999; (6) Strohmayer 
et al. 1996; (7) Mendez & van der Klis 1999; (8) Smith 
et al. 1997 (9) Wijnands k van der Klis 1997; (10) Boirin 
et al. 2000; (11) this work. 



